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Chapter 1 
Introduction 


The machine-sensible version of the AE-8 electron model environment was completed in 
December 1983. It has been sent to users on the model environment distribution list and is 
made available to new users by the National Space Science Data Center (NSSDC). AE-8 is the 
last in a series of terrestrial trapped radiation models that includes eight proton and eight 
electron versions. With the exception of AE-8. all of these models were documented in formal 
reports as well as being available In a machine-sensible form. The purpose of this report Is to 
complete the documentation, finally, for AE-8 so that users can understand its construction 
and see the comparison of the model with the new data used, as well as with the AE-4 model. 
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Chapter 2 

Overview of the Model 

The model is comprised of three parts: ( 1) an inner zone that covers the L (Mcllwain 
parameter) range 1.2 - 2.4, (2) an outer zone that covers the L range 3-11, and (3) a transition 
region for 2.4 < L < 3.0. For each part there are two versions: One depicts solar maximum 
(SOLMAX) conditions, and the other is for solar minimum (SOLMIN). An overview of each part 
follows. 


Inner Zone 


The inner zone parts were produced and documented earlier as models in the overall series: 
AE-5 Projected (AE-5P) for SOLMIN (Teague and Vette, 1974) and AE-6 for SOLMAX (Teague et 
al., 1976). Both of these were based on AE-5 (Teague and Vette, 1972) and the revision of the 
OV3-3 flux levels used in AE-5. AE-5, with an epoch of October 1967, was developed in a 
unidirectional differential energy flux form and converted to omnidirectional integral energy 
flux for its final output form. Analytical functions for the pitch angle dependence, including 
an atmospheric cutoff, and the eneigy spectrum were chosen based on the data. Three time 
effects were studied and modeled: magnetic substorms. Starfish decay, and solar cycle. For 
electrons below 0.7 MeV, the main temporal effect was a result of the solar cycle. Above this 
energy and at the L values above 2, substorm injections were large but infrequent. At epoch 
the Starfish residue was important for L < 1.7 mainly at intermediate energies. AE-5P was 
obtained from AE-5 by using the Starfish decay model of Teague and Stassinopoulos (1972) to 
remove the Starfish residue and by adjusting the solar cycle dependent fluxes of AE-5 to 
SOLMIN. The AE-6 model consisted of the removal of the Starfish residue from AE-5; no solar 
cycle changes were needed, since the epoch of AE-5 was one of SOLMAX. A later study of a 
number of new inner zone electron data sets indicated that no changes to these models were 
needed (Teague et al., 1979). 
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Overview 

Transition Region 


This region is the most difficult to model since it is affected by infrequent but large substorm 
injections that result in greatly varying energy spectra. Long-term (~ yearly) averages do not 
provide the quasi-steady picture found in the outer zone. Since the analysis of data to produce 
models differs between the Inner zone and outer zone, the decision was made to use a transition 
region as an interpolation region between the two models. In practice this means prescribing 
the flux at specified energies on the five L shells 2.5, 2.6, 2.7, 2.8, and 2.9. This was done as 
part of the original AE-5P and AE-6 models to stitch them together with the AE-4 outer zone 
model. Naturally, the values for AE-8 in this region are different from these earlier models. 
These differences are mainly for energies above 1 MeV, as shown in Chapter 4. Both SOLMIN 
and SOLMAX transition region versions exist, since the solar cycle effect is seen from L =1.4 
out to L = 5.5. 


Outer Zone 


AE-8 in this region was formulated by 

• Studying the data from Hovestadt’s Proton-Electron Detector experiment flown on the 
German Research Satellite A (GRS-A); Vampola's Energy Spectra, Fluxes, and Pitch Angle 
Distribution of Electrons experiment flown on OV1-19; and Vampola's Magnetic Electron 
Spectrometer experiment flown on OV3-3. 

• Using the analyzed results of Paulikas and Blake (1979) of the energetic electron data from 
their experiment on ATS 6 and Mcllwain's experiment on ATS 5. 

• Using the outer zone AE-4 model (Singley and Vette, 1972a and 1972b) to make 
modifications dictated by the data above to arrive at the outer zone portion of AE-8. 

The GRS-A satellite is commonly called Azur; this latter name will be used in this report. 

Using the Azur data, an empirical atmospheric cutoff was determined for AE-8 to replace the 
conservatively (pessimistically) chosen h max = 200 km for AE-4. The local time variation was 
extended down from L = 5 to L =3, also using the Azur data. Based on the ATS 6 local time 
variation obtained from Paulikas (1981), no solar cycle effect for this variation exists in AE-8, 
as opposed to AE-4. The energy spectrum above 2.5 MeV was enhanced at L = 6.6, using the 
geostationary data and then using Azur and OV1-19 data, the spectra above 2.5 MeV were 
altered from AE-4 throughout the outer zone. The flux variation along the field line has the 
same dependence as AE-4, and the statistical variation given by the standard deviation of the 
logarithm of the flux is also unchanged, except for the extension to higher energies. 
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Chapter 3 

Data Processing/Analysis 

The framework for the outer zone portion of AE-8 will be discussed now so that the processing 
and analysis of the new data can be understood. The details of the inner zone formalism will 
be reviewed in Chapter 5 (page 5-3 through 5-5). The formalism for AE-8 outer zone (OZ) is the 
same as that used for AE-4, where the omnidirectional flux above an energy (E ) is given as a 
function of the position along a field line ( b ) , the L shell ( L ) , the local time (<()), and the 
epoch ( T ) . The variable b is 

b = B/B 0 (3.1) 

where B c is the magnetic field strength at the geomagnetic equator based on an internal field 
only. B is the magnetic field strength of the internal field at the position in question; thus b is 
a coordinate along the field line that has the value 1 at the equator and °° at the dipole. The 
local time variation is mainly the effect of the magnetic field generated by current sources 
external to the Earth's surface. Although this choice of variables for the model is not the best 
for studying the details of particle dynamics in the outer zone, it is much more convenient for 
the user to apply the model to satellite orbits. 

The AE-80Z omnidirectional flux is written as 

J (>E, b, L, <{>, T) = N t (>E, L) <I> (>E, L, <{» G (b, L) (3.2) 

where Nx gives the equatorial integral spectrum, O gives the local time variation, and G 
provides the variation of flux along the field line. This simplification of three product 
functions is justified by the outer zone data. Specifically, the local time function has the form 

r w«|)-ii) i 

0(>E,L,<j» = K(>E,L) 10 L C(>E ' L)cos ^ J ( 3.3) 
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Data 

The peak In this function at <j> = 1 1 hours is an approximation that is good down to 0.2 MeV, 
and it was felt, as in AE-4, that more detail was not necessary. The function K is merely one of 
normalization so that 

2n 

Jtt = 1 (3 - 4) 

0 

The function C represents the amplitude of the local time variation as seen in the log of the 
flux. In AE-4 O was dependent on epoch; this is no longer true for AE-80Z as noted earlier and 
discussed further in the next section. 

The G function has the value of unity at the geomagnetic equator and is 

-mip fbc— q> +0.5 

Ibc-lJ 

G(b,L) = 

0 

The value for the exponent m is the same as AE-4, since this quantity can best be determined by 
data covering the low and mid latitudes, not high latitudes. The cutoff given by b c can only be 
obtained from high latitude satellites such as Azur. It is useful to note that the form of the b 
dependence given by G for an omnidirectional flux has an analytical counterpart for the 
unidirectional flux perpendicular to the magnetic field line given by (Roberts, 1965) 

— b^ — toK/jr . b < 

27l|i<0.5,m + l)lbc-lJ Vb c -1 ' ° 

9 p (b, L) - 

0 ; b > 

where f}(0.5, m+1 ) is the well known P function. Equation (3.6) was used in OV1- 19 data 
analysis. 

The processing of the Azur data will be described next. This will be followed by the OV1- 19 
processing, then that for OV3-3. Finally, the ATS 5 and 6 analyzed data will be discussed. 


( 3 . 6 ) 


; b < b c 
; b > b c 


( 3 . 5 ) 


Azur Proton-Electron Detector 


Azur was launched on November 10. 1969, into a 103' inclination orbit with perigee and apogee 
altitudes of 387 km and 3150 km. The Proton- Electron Detector experiment consisted of two 
cubical silicon (lithium-drifted) detectors, each covered by hemispherical shielding. There 
were two pulse height discriminators associated with each detector. This instrument provided 
four channels of information: (1) electrons greater than 1.5 MeV and protons greater than 20 
MeV, ( 2 ) protons 20-45 MeV, (3) electrons greater than - 4 MeV and protons greater than 40 MeV, 
and (4) protons greater than 72 MeV. This experiment provided data over the lifetime of the 
mission and covered the time period November 25, 1969 - March 15. 1970. The original data set 
supplied to NSSDC (ID# 69-097A-04A) contains count rates (CR) from the four channels 
described above along with all the position, time, and correlative data needed to analyze the 
data. More details on the instrument can be found in Achtermann et al. (1970). Ebel (1972) has 
given the details of the data processing of the raw data by the experimenters. 
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Data 


For L > 3, channels (1) and (3) above were measuring only electron fluxes. The conversion 
from CR to flux was provided by the experimenter. The calibration curves were not provided, so 
the usual threshold energy/ average efficiency studies outlined by Slngley and Vette (1972a) 
could not be done. However, comparison with other data showed that the threshold of 1.5 MeV 
was good while the higher energy threshold was better suited to 4.5 MeV for analysis purposes. 

The data were binned Into selected L shells of ± 0.05 width and tagged by the geographic local 
time <t> and the value of b. Plots were made of the flux versus b and. using equation (3.5) with 
m values from AE-4, a value for b c was determined. Having this, each CR at b = b i and L = 
Lj was projected to the equator by using the G function. This procedure is given by 

A = CR (bi, Lj) /G (bi, b^ j ) (3.7) 

and since G Is unity at the equator (b = 1 ) , one has 

A = CR (b = 1, Lj) (3.8) 

Then these data projected to the equator were binned In one-hour <|) cells to determine the local 
time distribution. The empirical distribution is given in Table 1 with a normalization the 
same as that given In equation (3.4). Any bins with no data were Ignored, and the integral (or 
sum In this case) was divided by the number of occupied bins. The least squares fit to the log of 
the flux provides values of C, which are also given In Table 1. The 95% confidence Interval 
half width for C, which is quite small, is also given there. The value of K was also determined 
in the fitting process, and, as expected, its value given in Table 1 is close to unity. In any 
satellite orbit except geostationary (or near geostationary), it requires a long time to sample all 
the local time bins. The coverage for Azur, which lived for about four months, can be seen in 
Table 1. A highly elliptical satellite such as OGO 1, 3, or 5 takes one year to sample all local 
times; geostationary satellites require only 24 hours. Consequently, there is "noise" 
introduced in the local time "signal" by enhancements (substorm Injections) in the outer zone. 
One expects that the local time variation should cease around L = 3. The values of C at 
energies above 1.5 MeV for various L values are plotted in Figure 1 for AE-4MAX and MIN, for 
Azur, and for ATS 6. The interval for ATS 6 data corresponds to SOLMIN conditions. The 
curve for AE-8 is shown; this applies to all solar cycle conditions. The AE-4MIN points are 
now believed to be elevated because of "noise" in the elliptical satellite data. AE-4MAX had the 
benefit of having the ATS 1 data of Paulikas and Blake incorporated. Below L = 6 the AE-8 
curve differs from the AE-4MAX model so as to conform to the Azur results. The fall of the 
curve at high L values reflects the fact that as the flux goes to zero (literally to J =1, which is 
about the value of the cosmic-ray proton background), there is no local time variation. 

After the AE-8 function was determined, it was used to convert the Azur data base to the 
local time of <|)m (L) • which is defined by the equation 

= 1 (3.9) 

and explicitly is 

12 K 

<t>m = 11 + — cos-1 (-“) (3.10) 

7t 

If a measurement was made at bi , L j , <j>i , then to convert this to (J^ j the following 
procedure was used. 

CR (bi f L j r ) = CR (bi, L j t <J>i ) /O ( L j , <j>i ) (3.11) 
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Data 

For each L value this converts the measurement to that which would be observed for the local 
time average at bj.. The CRs were then converted to flux using the geometrical conversion 
factor provided by the experimenter. Then the data were plotted with J versus b at various L 
values. The results for L = 3.0. 4.0, 5.0, and 6.0 are shown in Figures 2 through 5. In this region 
the data provide good coverage where the cutoff effects occur. The curves that are drawn require 
some explanation. These have the form of equation (3.5) with a multiplicative factor inserted 
such that equatorial value of J is the time average of the Azur fluxes projected there, as 
discussed earlier. The two curves in each figure differ only in the value of b c that was chosen. 
The larger value was chosen so that practically all the points have a b < b C 2 - The smaller 
value was chosen to represent a boundary for "undisturbed" conditions. This is best seen in 
Figure 2 where there seems to be a boundary of some small width with a relatively small 
number of points with larger b values. This curve was selected to pass through the center of 
the boundary, not its outer edge. As one looks at the higher L values, one sees a much higher 
percentage of points with b values above b c i. A vertical dotted line at b c i is included to help 
the reader see the points beyond this undisturbed boundary. The figures also illustrate why the 
Azur data are not good for determining the m value of equation (3.5). Satellites with 30' - 50' 
inclinations provide good coverage of the b region where cutoff effects are not present; all the 
available data of this type were analyzed in the construction of AE-4. 

The b C 2 cutoffs are the ones used in AE-8. The reason for this choice was to ensure that one is 
very unlikely to experience fluxes beyond the AE-8 cutoff. Feedback from the users have 
indicated this is preferable to them. The b c i cutoff is more useful in studying the particle 
population. This is illustrated in Figure 6 where the cutoff values are plotted versus L. The 
inner zone values are taken from the AE-5 document (Teague and Vette, 1972). Notice that the 
outer zone b c i is virtually an extension of this inner zone boundary with nearly the same L 
power law dependence. In the outer zone b C 2 is about a factor of 1.4 greater than b c i. After 
passing through the transition region, there is no distinction between these two boundaries in 
the inner zone. 

The greater than 4.5 MeV channel was processed in the same manner as described earlier in 
this section for the greater than 1.5 MeV data. The data from this channel showed the same 
atmospheric cutoff (b C 2 ) as the greater than 1.5 MeV data. The dynamic range for the greater 
than 4.5 MeV channel was not as large as that for the greater than 1.5 MeV, being about a factor 
of 50, and the intensity threshold corresponded to a flux of about 20 electrons/cm 2 -s. Because 
of these limitations the data did not reveal a consistent local time distribution. So the C (> 1.5 
MeV) was scaled using the AE-4 model energy dependence to obtain the local time distribution. 

Returning to Figures 2 through 5. if one imagines moving the curve associated with b C 2 so that 
it passes through each plotted point in chronological order, then the projected equatorial flux 
value as a function of time is obtained. The daily averages of these projected equatorial fluxes 
from the two Azur channels are shown in Figures 7 and 8 for the L values corresponding to the 
outer zone peak. The daily sample size varied from 5 to 18 and averaged about 10. The 
standard deviations of the logarithm of the flux (a) for these plots are 0.51 and 0. 12, 
respectively. The statistical function for AE-8, which uses a is discussed in Chapter 5 (pages 5- 
5 through 5-7). It was felt that the new data used in AE-8 could not add much to this statistical 
function over that provided by the data used in AE-4. Both a values above from Azur are low 
compared to AE-4. Data sets that have a duration of less than six months often show lower (or 
higher) as than longer lived ones. The use of daily averages also reduces the a value. It is 
concluded that Azur statistical behavior is similar to other outer zone data sets. 

Thus, the Azur data have supplied the omnidirectional flux for E > 1.5 and > 4.5 MeV for L 
values from 3.0 to 6.6. It is believed that the greater than 4.5 MeV data suffer from a 
bremsstrahlung background that becomes important for L above 5.5. Another contributing 
factor might be the intensity threshold for this channel; when flux levels fall more often below 
this threshold, the average value is distorted upwards. This effect would also reduce the value 
of the a just discussed. 


3-4 



Data 


OV1-19 Electron Spectrometers 


OVl-19 was launched on March 18, 1969, Into a 105° inclination orbit with perigee and apogee 
altitudes of 466 km and 5764 km. Thus, it was in a similar but higher orbit than that of Azur. 
Vampola had an experiment on board that was comprised of two magnetic electron 
spectrometers. These instruments provided data from March 18, 1969, to the end of spacecraft 
life in March 1970. The data set deposited in NSSDC (ID# 69-025C-05A) terminates on January 
25, 1970; consequently, there is a data overlap period with Azur of about 61 days, neglecting 
data gaps. 

The initial results from this experiment were used to claim that the AE-4 model was in gross 
error for electrons whose energies are above 1.5 MeV. Unfortunately, it was many years after 
these distorted (but not completely false) claims before the resolution of this matter was 
obtained. This came about through the analysis described here, the result of which was the 
utilization of much of the data from this experiment in AE-8. However, the channels above 2 
MeV are rather useless in the outer zone as independent data, since a background correction is 
absolutely necessary. 

The description of the instrument was given by Vampola (1971) and will be summarized here 
since an understanding of it is necessary to follow some of the analysis. The spacecraft was 
spin stabilized at a rate of 8.4 rpm. Both spectrometers had their look direction perpendicular 
to the spin axis. A three-axis magnetometer was also on board. The low energy spectrometer 
(LES) had eight detectors for measuring electrons and one detector for background, all in an 
analyzing field of 470 gauss. The background detector had two channels, one for measuring 
protons greater than 55 MeV and one for measuring electron bremsstrahlung (BREM). The 
BREM discriminator was set to the nominal energy loss level used by the eight detectors. The 
high energy spectrometer (HES) contained 16 detectors plus a background detector with two 
channels, one for protons greater than 130 MeV and one for electron BREM; this latter 
discriminator was set in the same manner as that for the LES. All detectors were lithium- 
drifted silicon 1 mm thick with an unshielded sensitive area of 0.95 cm * 1.37 cm. The energy 
levels and other information about this instrument are given in Table 2. Both spectrometers 
were the 180° focusing type with internal baffling around the edges of the uniform field and on 
the pole plates to provide a disk-loaded collimator. A tungsten collimator was used at the 
entrance of the magnet. Magnet pole pieces were made of Index V, and an iron yoke enclosed 
the whole assembly. From the threshold energy of protons arriving at the background 
detectors, one can infer that the HES had about five times as much shielding as the LES. 

Data on magnetic tape (NSSDC ID# 69-025C-05A) contained the CRs for each detector along 
with all the ancillary data needed for analysis. A subset of these data was selected for AE-8 
analysis; only data with local pitch angles within 3° of the perpendicular to the measured 
magnetic field as determined by the on-board magnetometer were studied. Thus, one is dealing 
with the unidirectional flux perpendicular to the field line, denoted here by j p . The CR meters 
were sampled every second: however, only eight detectors plus two background channels were 
read out at one time. LES was read and then two reads were required for the HES. 

The big mystery for this experiment was the background. The proton channels had such small 
counting rates that they could be ignored in the outer zone. However, the HES BREM channel 
counted much too high to be background for the HES detectors. In Table 3 the ratio (CR - 
BREMJ/CR, averaged over the whole interval March 18, 1969, to January 25, 1970, is given. 
Negative ratios mean that the background exceeded the raw electron CRs, and from these 
numbers one can rule out the use of all channels above 1.7 MeV and the 53 keV channel over the 
range of L values shown, unless appropriate backgrounds can be obtained for subtraction. The 
ratio of the two BREM channels has also been studied without shedding any light on the 
problem. Ignoring the BREM channel, as Vampola did, one arrives at flux values that are much 
too high (factor of 10-100) compared to other data, as shown in Chapter 4. 
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Data 


The ATS 6 data of Paulikas and Blake (1979) showed clearly in the geostationary region for 
energies above 2.5 MeV that AE-4 was too low and therefore needed revision in the high energy 
range. Since the OV1 - 19 data was the only data set that covered the energy range between 2 and 
5 MeV at other outer zone L values, it was imperative to try to use these data to revise the 
outer zone . 

The following technique was developed. In Figures 9 through 1 1 the OV1- 19 data are shown 
projected to the equator using the same techniques applied to Azur data except that equation 
(3.6) was used instead of (3.5) . The open solid squares, labeled uncorrected, are the raw CRs 
converted to flux using the conversion factors of Table 2. Notice that the spectrum is flat above 

about 3 MeV for all the L values. This raw directional flux, j r , consists of the true electron 

flux, j p , plus a background flux, j^, resulting from BREM; i.e., 

jr = jp + jb (3 . 12) 

By taking the last three channels of the HES with the background (supposedly) removed and 
performing the following operation 

S = j p (5091 keV) [5243-4942] + j p ( 4 7 83 keV) (4942-4633) + 

j p (4476 keV) [4633-4500] (3.13) 

the total electrons/ cm 2 -ster-s in the energy interval 4500 to 5243 keV results. Now this can be 
converted to an omnidirectional flux at the equator In electrons/ cm 2 -s by multiplying by the 
ratio of equation (3.5) to (3.6) 

f = G/g p = 2Jtp(0.5, m+1 ) (l-l/b c ) 0 - 5 (3.14) 

The product f S obtained from equations (3. 13) and (3. 14) is denoted here by J e 3 - Then 

Je3 s J (> 4500 keV) = Ja(> 4500 keV) (3.15) 

where J A is the Azur data. This is approximately true since J (> 5243 keV) in the outer zone is 
small relative to J (> 4500 keV). If j r is substituted for j p in equation (3.13), then J r 3 
replaces J e 3 in equation (3.15). J r 3 is larger than J A , the amount depending on the 
background contamination. Taking the ratio 


d “ Jjj < 3 ' 16 > 

the fraction of the uncorrected flux representing the real electron flux results. The empirical 
results for d are shown in Table 4. The data from both the Azur and OV1- 19 experiments were 
averaged over the data overlap period November 25, 1969 - January 25. 1970, to obtain this 
ratio. It is quite clear that in this energy range the OV 1-19 experiment was highly contaminated 
with BREM from lower energy electrons. This quantity d has the same meaning as the ratio 
given in Table 3 except that the background is determined by using the Azur > 4500 keV data. 

The BREM channels had some 25 to 90 times more counts than expected. The most likely cause 
of this was that the HES BREM discriminator level in orbit was much lower than those for the 
HE1-16 detectors. In the case of the LES, the problem may have been that the LEI and LE2 
discriminator levels were set high relative to the LES BREM discriminator. The considerably 
lower conversion factor in Table 2 suggests this was the case and is so noted in the remarks 
column. This would mean that the background for these two channels would be less than that 
Indicated by the LES BREM channel. Noting that 
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• The six highest channels of the HES showed nearly the same unconnected flux for 
the L Interval 3.0 to 6.6, 

• Most of the uncorrected flux was produced by BREM, 

• The HES BREM channel had much higher CRs than many of the HES detectors and 
could not be used to properly subtract the background, 

the following procedure was used to determine the average HES background. This is expressed 
as a directional background flux, j^. This was obtained by averaging the six highest HES 
channels and applying the factor ( 1 -d ) 

HE1 6 

jb = £jr(i) (3.17) 

i=HEl 1 

Since Figures 9 through 1 1 show that these j r s are essentially independent of energy, it is 
believed that jb applies to all of the HES detectors. Then using equation (3. 12) the value of j p 
can be obtained from the uncorrected flux. In the few instances where this procedure resulted 
in a negative j p , then the equation 

j p (k) = dj r (k) (3.18) 

was used. These results are plotted as the open squares, labeled corrected, in Figures 9 through 
11. It is apparent that the corrected and uncorrected data above 3 MeV differ by more than a 
factor of 10. The results also Indicate that no correction is required for the LES detectors. 

The first corrected Integral spectrum for OV1-19 was then constructed by 

17 

J V <> Ei) = f£j p (k)W(k) + J A (> 4500 keV) 
k=i 

+ 175df j r (4476 keV) (3.19) 

where W (k) Is the energy channel width for the k th useful channel, and the factor 175 (keV) is 
the width from 4325 to 4500 keV. The useful channel Index and the Ws are identified in Table 5. 
In the final results for OV1-19, J AE _ 8 (> 4500 keV) was used in equation (3.19) in place of J A (> 
4500 keV). 

The numerical results of all these procedures of this section are given in Tables 6A through 6D. 
The comparisons of the final corrected data with other data are presented in Chapter 4. 


OV3-3 Electron Spectrometer 


The OV3-3 satellite was launched into a 81.4° inclination orbit with perigee and apogee 
altitudes of 360 km and 4492 km. Vampola had a magnetic electron spectrometer on board 
very similar to those on OV1-19. This one used nine detectors for measuring electrons, and one 
of them had a discriminator channel for measuring penetrating protons greater than 105 MeV. 
Since BREM was found not to be a problem in calibration, no BREM background channel was 
provided. The detectors, collimation, and baffling were nearly identical to those used later in 
OV1-19. Each detector had an independent set of electronics consisting of an amplifier, two- 
level pulse height analyzer and log(CR) to analog converters covering the CR ranges 1-10 3 and 
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300-2 x 10 5 . More details have been given by Vampola (1969). The channel energy widths and 
geometric factors are given In Table 7; also Included are the threshold energies for Integral 
spectrum calculations that are employed In the data processing done for AE-8. The geometric 
factors are four times greater than those published by Vampola and are based on his revision 
(Vampola. 1972, as given In Teague and Vette, 1974). 

The data were supplied as flux values versus time and Included the pitch angle as Inferred from 
a model geomagnetic field along with other parameters needed for analysis. The subset of this 
data used In the AE-8 analysis consisted of the data for pitch angles close to 90° for L < 10. The 
NSSDC data set (ID# 66-070A-05C) covers the time period August 4, 1966, to September 9, 1967. 

The main purpose In processing the OV3-3 data in the outer zone was to provide a comparison 
with the OV1-19 data In the energy ranges that overlapped. The OV3-3 data had been used In 
producing AE-5, and when Vampola changed the geometric factors of the Instrument, it was 
necessary to produce AE-6 to replace AE-5. The processing for OV3-3 In the outer zone was done 
in an identical manner to that used for OV1- 19 with one exception. The Integral spectrum J 
(> 2.445 MeV) was estimated by extrapolating the OV3-3 differential spectrum to higher 
energies. This procedure was not possible for OV1-19 because of the background problem. The 
results are shown and discussed in Chapter 4. The agreement between the OV3-3 and the 
corrected OV1- 19 data Is good at certain L values. Time variations can explain the 
disagreements at other L values. 


ATS 5 and 6 Omnidirectional Spectrometers 


On ATS 5 Mcllwain flew an Omnidirectional High-Energy Particle Detector experiment 
comprised of three plastic scintillator detectors with a 2rc field of view and hemispherical 
shielding. The Instrument measured electrons In the range 0.5 to 5 MeV in 12 channels and 
protons greater than 12, 16, and 24 MeV. Paulikas flew an Omnidirectional Spectrometer 
experiment on ATS 6 that consisted of four solid state instruments. One was a two-element 
telescope with a 30° full cone angle, and the other three were omnidirectional Instruments 
similar to the Azur detectors described earlier. Particles and energy ranges (In MeV) measured 
were 


Electrons: 0.14-0.60, > 0.7, > 1.55, and > 3.9 

Protons: 2.3-5.3, 3.4-5.3, 12-26, 20-52, and 40-90 

Alphas: 9.4-21.2, 13.4-21.2, and 46-100 

The lowest energy electron channel and the two lowest proton and alpha channels were 
directional (from the telescope); all the rest were omnidirectional. Only the electron channels 
of these two experiments are of interest here. 

Paulikas and Blake (1979) presented the data from Mcllwain's ATS 5 experiment and their ATS 
6 data in the following manner, as shown In Figure 12. They produced from Mcllwain's data 
the omnidirectional fluxes to match their four channels on ATS 6. The ATS 5 data cover the 
approximate period 1969.7-1972.3 while the ATS 6 data cover 1974.5-1978.3. The ATS 5 data 
were normalized to the ATS 6 data during some simultaneous coverage in mid- 1974. The 
Zurich monthly sunspot number is also plotted in Figure 12 and shows that the combined data 
cover both SOLMIN and SOLMAX conditions. The ATS 1 data of Paulikas used In AE-4 also 
appear In Figure 12. All of the electron data are plotted as 27-day running averages, and yearly 
averages are shown by the horizontal bars extending across each year or partial year. These 
yearly or partlcil year averages provide the pertinent parameters for AE-8. It can be seen In 
Figure 12 that there are no trends In these averages associated with the solar cycle. The 
variations are less than a factor of two except for the highest energy channel where variations 
as large as a factor of 5 are seen; these are more stochastic In nature and possibly related to 
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high speed solar wind streams according to Paulikas and Blake (1979). The three satellites 
were at slightly different coordinates, since they were operating at different longitudes. ATS 1 
was right on the geomagnetic equator at L = 6.6 while ATS 5 and ATS 6 were slightly off the 
equator at L - 6.7 and 6.8, respectively. Paulikas and Blake ignored these differences since they 
were unimportant to their study. 
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Chapter 4 

Comparison of the AIDS Model 
with Data 


This comparison is limited to the data discussed in Chapter 3 and to the AE-4 model. The 
documentation for AE-4 and AE-5 gives the comparison of many data sets to the models, and 
the work of Teague et al. (1979) shows comparison of other data sets with inner zone models. 
Therefore, this section is confined to outer zone electrons, whose phenomenology has been 
discussed for modeling purposes by Singley and Vette (1972b) and by Chan et al. (1979). These 
comparisons are shown in Figures 13 through 39. In all of these AE-8MAX (MIN) is represented 
by the large open squares (circles) and the other entries (except for Figure 33) by smaller 
symbols so that the reader can clearly see the data to which the new model conforms. 

Figures 13 through 17 cover the transition region. Each is broken into two graphs: (a) SOLMAX 
and (b) SOLMIN. For SOLMAX, AE-8 and AE-6 agree up to energies of 0.75 MeV. Above this 
energy AE-8 becomes larger and extends to higher energies. For SOLMIN the agreement 
between AE-8 and AE-5P is good up to 3.5 MeV at which point AE-8 grows larger and extends to 
higher energies. In Figures 17a and 17b AE-4MAX and MIN replace AE-6 and AE-5P, 
respectively. Here for SOLMAX the agreement extends up to energies of 3 MeV, then AE-4 falls 
rapidly below AE-8. For SOLMIN the departure of AE-4 from AE-8 begins at 0.75 MeV, with AE- 
4 falling rapidly below AE-8: again AE-8 extends to higher energies. 

The next set of figures that will be discussed involves Figures 18 through 32. In these the 
comparisons of AE-8MAX with AE-4MAX and data from Azur, OV1-19 unconnected (19U), OV1- 
19 corrected (19C), and OV3-3 are presented. OV3-3 is presented only for L values of 3.0, 3.5, 

4.0, 4.5, 5.5, and 6.0; thus, these data are only found in Figures 19, 20, 22, 23, 25, 26, 28, and 29. 

Some remarks relative to this set of plots are appropriate. First, the focus will be on Figures 18 
through 21, which cover the L region 3.0 through 3.5. The electron fluxes above 0.5 MeV in this 
region are very volatile, even when one considers yearly averages. This is the same volatility 
as that demonstrated in Figure 12 for the greater than 3.9 MeV fluxes. If one concentrates on 
the energy region 0.75 through 2 MeV, one finds at L = 3.0 that 19C is low, Azur Is in the 
middle, and OV3-3 is high with the ratio OV3-3/ 19C - 20. At L = 3.5 Azur is low, OV3-3 is 
middle, and 19C is high with the ratio 19C/Azur - 6. The reason for this behavior has been 
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discussed previously (Chan et al., 1979) and Is related to Injection events that penetrate deeply 
Into the radiation belts. During the lifetimes of OV3-3 and OV1-19, there were two major 
injection events that penetrated to L = 3.5, but only one of those four events penetrated to 
L = 3.0. This latter event occurred during the OV3-3 lifetime. During the Azur lifetime only 
one of the OV1- 19 events was present. Consequently, the observed behavior seen in Figures 18 
through 2 1 Is understood. The position of AE-8MAX is consistent with AE-4MAX and lies near 
the geometric mean of OV3-3 and 19C at L = 3.0. The high eneigy extension above 2.5 MeV for 
AE-8MAX follows the Azur and 19C results. It is noted that there Is forced agreement between 
19C and AE-8MAX at 4.5 MeV, as explained in Chapter 3 (page 3-5 and subsequent pages). Below 
0.5 MeV there is relatively good agreement. The extreme of the 19U data above 2 MeV is obvious 
throughout this L region. 

A few remarks relative to errors are appropriate at this time. The general error assigned to the 
models in this series has been stated as about a factor of 2. For AE-5 a more detailed statement 
was made. What Is meant by the factor of 2 error is that for a flux of J m the reasonable upper 
limit Is J u and the reasonable lower limit is Ji where 



Thus, the ratio = 4. Now the symbol for AE-8 in Figures 18 and 19 would have to have a 

vertical extent about three times longer (linear not the log scale) to represent this error. The 
actual dimension represents a factor of 1.25. Consequently, the symbol can be used to roughly 
gauge how well the data and the model agree. Therefore, in the worst situation at L = 3 the 
error is a factor of 4.5. Elsewhere in this region the factor of 2 is appropriate and, in some 
cases. Is pessimistic. Only data lying outside the factor of 2 error will be noted for the 
remainder of the discussion of these figures. 


In Figures 22 through 32 all is good save for a few spurious points until the Azur greater than 
4.5 MeV data gets too high by L = 4.6 and continues this trend at the higher L values. The 
disparity of the 19U data is obvious throughout this region. It should be noted in Figure 32 that 
AE-8MAX and AE-4MAX are replaced by AE-8 and AE-4 since the solar cycle effect is no longer 
evident and the MAX and MIN models are the same at higher L values. 

For L = 6.6, the nominal geostationary L shell, instead of plotting AE-4, all of the satellite 
data used in constructing both AE-4 and AE-8 are shown. This requires two figures because of 
the large number of satellites involved. In Figure 33a the ATS 5/6 data are plotted for both the 
highest (U) and lowest (L) averages shown in Figure 12. It should be noted that the 4 MeV model 
value conforms to ATS 5/6 U value and so is a pessimistic value. At least the neighboring 19C 
points are in accordance. The Azur point at 4.5 MeV is clearly too high, as discussed in Chapter 
3 (page 3-2 and subsequent pages). The OGO 1 point around 2 MeV is high, and the ATS 5/6 point 
at 0.7 MeV is at the factor of 2 limit. The 19C points between 0.8, and 3 MeV are also errant. In 
Figure 33b the high eneigy points for Explorers 6 and 12 are low, but as ATS 5/6 data at high 
energies have shown, fluctuations in long-term averages can be expected in this energy-L 
region. The high energy 19U points reach their ultimate departure from reality. For the latest 
ideas on time variations of high energy electrons in this region, consult Baker et al. (1986, 

1989, and 1990) and Nagai ( 1 988). 

For the higher L values 7 through 1 1 there are no new data to compare with AE-8. 
Consequently, in Figures 34 through 39, the comparison of AE-8 and AE-4 is given. It is 
apparent that AE-8 possesses a high energy tail except for L =11. 
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Chapter 5 

Presentation of the Model 

This model has several different forms of presentation. The most important for use is the 
machine-sensible or numerical form, which will be given first. Based on the analyses 
conducted to construct the model, there were certain analytical descriptions not incorporated 
in the numerical model that might be useful to some recipients; these are summarized next. 
Graphics provide a good way of conveying the model. However, some of the standard forms 
used in previous model documents do not appear, such as the b-L and R-A flux contour plots 
and the J versus b. The radial profile equatorial flux plots are included to illustrate the solar 
cycle effect. Finally, the results of fluences obtained in the standard set of circular orbits are 
given in tabular form. 


Numerical 

The function 

Z = log! oil 


where 


2n 

iI(>E,b,L,T) = J J(>E,b,L,<t>,T)d<]> 

0 

is given at a fixed set of instances of the variables. This presentation consists of both the 
logical and physical organization, since the physical format may vary depending on the 
storage medium on which the model is conveyed. 
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Logical Organization 


The logical organization for the numerical AE-8 was dictated by the desire to store the 
essential data in a minimum amount of storage space. This goes back to the early 1970s when 
the memory size of most computers was well under a megabyte. AE-4 and -5 were Just coming 
out and were larger than earlier models. Consequently, a BLOCK DATA statement approach 
was taken using binary coded decimal (BCD) representation with six characters per storage 
cell. Also, the integer format made for easier interchange; this required the scaling of the data 
by differing scale factors. 

The logical record chosen has been explained in Teague et al. (1972), but it will be described in 
detail here, since some updating is necessary. This logical record is variable in length with the 
hret field giving the number of fields In the record. The second field is the scaled energy value. 
The rest of the fields are organized first by L value and then by a variable related to b. In 
addition, the variable Z of equation (5. 1) is given at the equator and denoted Zi here The 
first field in the L substructure is the number of cells in the substructure followed by the 
scaled L value and then the field for scaled Z v The remainder of the cells in this substructure 
contain scaled values of 


5(b-l) = 8 ( x ) = 5b = 5x 

where the change is defined by a fixed Z decrement, ZDEL, which for AE-8 has the value of 
0.25. Therefore, on the L shell to reach 

z 2 = Zi — ZDEL 

the b value must change from 1 to 

bi = 1 + dbi 


Thus, the vaiue of 8bi is variable, but given Zi the value of Z i+ i is known from the cell 
index of Obi. This pattern continues for each of the L values in the AE-8 grid, which can total 
up to 45 depending on the energy. (The variable x is convenient for the method of 
interpolation described in Appendix A.) This logical record then repeats for the other 17 
energies (18 total) of the AE-8 grid. As shown in Table 8 this scheme utilizes 13,548 cells for 
AE-8MAX and 13, 168 for AE-8MIN. What are not contained in the basic model map are the 
scale factors for E, L, Z, and db as well as the total number of cells, ZDEL, name, and epoch, 
being redundant to the name. These are given in a header record and are included in 
Table 8. The parameter ZDEL does not appear in the BLOCK DATA statement. 


Physical Representation 


Magnetic Tape. The data set NSSDC ID# MT- 2AB contains AE-8MAX with a BCD 
representation and a physical record size of 80 bytes. The physical record is a card image- the 
header record and first part of the model map are given in Table 9. There are 1692 physical 
records in this data set. The first 28 physical records are header records, and there is one 
trailer record containing the END statement. The format of the BLOCK DATA is given in Table 
9. The model is broken into 84 (sub) maps with the first 83 being of dimension 163 and the last 
19. By decoding MAPI, one can see that the first two L values are 0 and 1. 1 with Zi = 0 . These 
are guard values at the inner edge of the map. At certain energies the L value of 1 1 is replaced 
by one of the grid values in Table 9 that is just prior to the L grid value for which Z has a non- 
zero value. The same two "empty" L sets occur at L values of 12.0 and 15.6 at the outer 
boundary, with the L value of 12 being replaced by a lower value for certain energies. The 
details of the L grid can be seen more clearly in Chapter 5 (page 5-7) where Table 13 is 
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discussed. AE-8MIN differs only in its total length being 1646 physical records, containing 79 
(sub) maps of dimension 163, one of 159, and the last of 132. The data set ID# is MT-2AA. As 
noted above, the value of ZDEL is missing from the BLOCK DATA statement. This parameter 
is contained in the software program MODEL sent with this data set and used to decode the map 
to produce flux values. 

Microfiche. These two data sets (AE-8MAX — #MT-52A and AE-8MIN — #MT-52B) provide 
integral flux values from L = 1.2 to 11.0 at 0.1 Intervals for B values (magnetic field 
intensities) from the equator to the flux cutoff limit of 1 electron/cm 2 -s in 100 gamma steps 
starting at the first multiple of 100 gammas. The energy levels covered are slightly different 
than the MAP energies being in MeV: 0.04, 0. 10, then the range 0.25-4.00 in 0.25 steps, and 
finally 4.00-7.00 in 0.50 steps. 

Floppy Disk. In 1988 Bilitza completed the revision of Program MODEL (MODEL87) and wrote 
an interactive driver program RADBELT to handle the models in the BLOCK DATA format 
(Bilitza. 1990). Included as part of this package is AE-8. The whole package has a data set ID 
number of PT-1 IB and is available on two PC floppy disks. The parameter ZDEL is included in 
the BLOCK DATA statement in this form. 

Since the main element of decoding the BLOCK DATA format is contained in the FORTRAN 
function TRARA2, this is described in Appendix A with some additional comments added to 
the source listing to make it easier to follow the logic. 


Analytical Forms 

There are several analytical expressions that have been used in constructing the model, some 
of which are not incorporated in the numerical model. An example is the local time function 
given by equation (3.3); this function also contains tabular functions. These analytical/ 
tabular forms will be summarized in this section. 


Magnetic Cutoffs Resulting from the Atmosphere 

For AE-8 this takes the form 


0. 6572L 3 - 452 

f 

1.2 

< L < 2.4 


b C 2 = 0.196L 4 ' 878 

r 

2.4 

< L < 3.0 

(5.4) 

1 .4567L 3 ' 050 

f 

3.0 

< L 



Results from equation (5.4) were incorporated into the numerical model. A more physical 
cutoff was also obtained but not incorporated. This was given by 


0. 6572L 


3.452 


bci = 


1 .0523L 


3.050 


1.2 < L < 3.23 


3.23 < L 


(5.5) 
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Local Time Function 


The function C ( > E , L ) is needed to use equation (3.3); the function K (> E , L ) can be 
derived once C Is known. C Is given as a tabular function in Table 10. As noted in Chapter 3 
(page 3-3), C has no dependence on epoch. 


Flux Variation Along the Field Line 


In the Inner zone this was given as an equatorial pitch angle distribution, differential in 
energy, in the form 


j(E,a 0 ,L,T) = a(a 0 ,L,T)Eexp(-E/X(a 0 ,L,T) (5.6) 


where 


9 (tto/ L, T) 


a r (L,T) 


sinP (q 0 -q c ) 
sinP (CCe-Ctc) 


Otc £ (Xq < (Xg 


(5.7) 


a r (L,T) ; a e < a G < 90° 

Here the cutoff Is expressed by 


Otc 



The function X of equation (5.7) is given by 


(5.8) 


o 1 n H/y 

x r (L,T) . n ; a c < a c < a e 

sin n a e — 

X ( a Q , L , T ) = (5.9) 


x r (L,T) ; a e < a Q < 90° 

Now the functions a r , p, O^, X r , and n were given as tabular functions for AE-5. However, In 
constructing AE-5P and AE-6 these functions were not tabulated. The only tabular function 
remaining for magnetic field line variation is in the outer zone and Is the exponent m of 
equation (3.5); this Is given as a tabular function in Table 11. 

It may be Instructive to compare in more detail the differences between equations (3.5) and 
(5.6). Since there is no local time dependence in the inner zone, the O function does not appear 
In equation (5.6). To show the comparison better, integrate equation (5.6) over energy to obtain 
the Integral spectrum 


Diz (- > (Xg, L ( T) — aX^(^ + l)exp(— ^) 


(5.10) 
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Then by substituting 

b = sin -2 a 0 (5.11) 

the perpendicular flux 



\ 

f- E (rY /2 1 

V^e ) 

exp 

J 

x r J 


(5.12) 


results. It is not possible to convert this analytically into an omnidirectional flux. 
Consequently, the best that can be done is to convert equation (3.2) into its j p form using 
equation (3.6) to obtain 


Dpoz 


27t|3(0 .5,m+l) l b c 


'be-b K, / b c 
bc-lj V b c -l 


N t (> E,L) 


(5.13) 


In the case where X r « E and keeping only terms that convey functional dependence 


- (p-n) /2(. A / b~ L 
Dpiz ~ b V 1_ \ bc-l 

-m m 

j poz ~ b (b c -b) N T (> E,L) 


Eexp 


fWfr ) 

v be / 


(5.14) 

(5.15) 


A comparison of these two equations shows that the inner zone spectral function depends on b 
as well as E; both spectral functions also depend on L and T. All exponents depend on L and 
T. The main field line variation near the equator expressed by the first factor in each equation 
is the same; however, the behavior near cutoff is of a different form. In addition, the 
exponential in equation (5. 14) has a profound effect on the b variation. It is possible to use the 
form of the b variation given in equation (5. 13) or (5. 15) to fit the one given in equation (5.12) 
or (5. 14) at a given energy by adjusting the value of m, but then m is strongly E dependent. 
Since AE-4 and AE-5 were developed somewhat independently, a common form for the b 
variation was not enforced. 


This discussion shows where the character of the inner and outer zone portions of the model 
are different and why a transition region is necessary. 


Statistical Function 

The statistical model associated with AE-8 is the same as that for AE-4 but extended to higher 
energies. The time variations in the outer zone make the flux a random variable. All random 
variables in this section will be denoted with bold face symbols. The empirical cumulative 
distribution function is obtained by 
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P(J > J s ) 


Time that J > J 3 
Total time of sample 


(5.16) 


where J is the random variable representing the local time averaged flux, which is still 
subject to large excursions resulting from injections. To a reasonable approximation this 
function P can be represented by 


P (Z > Z s ) 


oo 



f (Z- H) 2 

J exp (- 2<j2 dZ) 

Zs 


(5.17) 


which is the well known normal distribution, and Z is defined by equation (5.1). When Z s = 
JO., then P =0.5. ForAE-8, O (> E,L) is a tabular function given in Table 12. This function 
follows the AE-4 function and is extended in energy. As the flux tends to zero, o also tends to 
zero. The value of pi is obtained using o and the long-term average value of the flux, as will 
now be shown. The probability density function q is 


q(Z) 


1 ( Z— Hr) 2 

= exp (- ) 

271 2c 


(5.18) 


It is convenient to change variables from Z to Y to obtain a random variable with mean of 
zero and standard deviation of 1 . With 


(z-n) 

Y = 

a 

and J = lOM-exp (koY) 

where k = In (10) = 2.30259 
one can compute the average value of J 
oo 


J* 


1 

V 271 

J T 

_oo 


Y 2 

Jexp ( - 2 dY) 


= lOM-exp 


2 r\2 


k z o 


oo 

r 


j 

-OO 




27t 


(Y+km 2 

exp ( - dY) 


(5.19) 


(5.20) 

(5.21) 


(5.22) 
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Since the integral is unity and the exponential can be written as a power of 10, the result is 

kG 2 

logJ* = H + ~Y~ (5.23) 

With AE-4 it was noticed that some users considered the relationship to be 
logJ* = ji 

which would result from treating the 50% probability of occurrence as being the average flux 
value. Since G is typically equal to 0.5, this produces an error of about a factor of 2. 

Although the statistical model for AE-8 is not directly incorporated into the machine-sensible 
model, it can be done using Table 12 along with some suitable coding. This has been done at the 
European Space Agency (ESA) European Space Technology Center (ESTEC), for example. 


Graphical Presentation 


In Figures 13 through 39 the integral spectral plots of equatorial omnidirectional flux versus 
energy were used to demonstrate the comparison among AE-8, AE-4, and various data. The 
equatorial omnidirectional flux versus the L value is now used to convey some features of the 
model and to illustrate the solar cycle effect in some detail. The data base for these plots was 
extracted from the BLOCK DATA statement and is printed as Table 13A (for AE-8MAX) and 
Table 13B (for AE-8MIN). An examination of these tables shows the position of the L values 
where Z\ is equal to zero. Besides the special values of L (0, 1. 1, 12, and 15.6) discussed earlier 
(page 5 - 2), the values 2.69, 2.79, and 2.89 also appear to handle the region near where the flux 
minima occur in the slot region. 

In Figure 40 the equatorial flux profile for all of the grid energies are shown for AE-8MAX. No 
electrons greater than about 4.5 MeV appear in the inner zone, while the outer zone contains 
electrons up to 7 MeV. When the energies are above about 0.75 MeV, the inner zone peak occurs 
at an L value of 1.50, while at lower energies this peak is at 1.80. For energies above 2.0 MeV 
the slot is very deep. In the outer zone the flux peak moves inward from about an L of 4.8 as the 
energy increases to about an L of 3. 1 for the highest energies. 

Figures 41 through 45 show the equatorial profiles for AE-8MAX and AE-8MIN for the same 
energies, so that the solar cycle effect is clearly demonstrated. Large open symbols have been 
chosen for AE-8MAX, while small solid symbols have been used for AE-8MIN. With this choice 
it is easy to see when the fluxes are essentially the same value. From these figures it can be 
ascertained 

• For energies below 0.75 MeV the inner zone fluxes are higher at SOLMAX. 

• For energies above 1.0 MeV the inner zone fluxes are independent of the solar cycle. 

• For energies above 6.0 MeV the outer zone fluxes are independent of the solar cycle. 

• For energies below 6.0 MeV there is a region starting at the outer edge of the slot and 
ending in the L range between 3.0 and 6.0 depending on energy where the SOLMIN 
fluxes are less than SOLMAX fluxes. The ratio of these fluxes ranges up to 4 at the lowest 
energies. 
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Orbital Integration Results 


There Is a standard set of circular orbits that are run with the models to demonstrate the 
amount of fluence In electrons/cm 2 -day accumulated In various energy windows and above 
certain threshold energies. The results are shown In Table 14 for AE-8MAX. The sampling 
Interval and the accumulation time vary with orbit altitude to take Into account certain 
factors and provide a reasonable average value. Thus, at low altitudes the sampling time Is 30 
seconds because the satellite moves rapidly through the regions where flux gradients are large. 
The accumulation length is 36 hours to average out the South Atlantic Anomaly effects. These 
parameters are shown at the top of each page of the table. For orbits that traverse the Inner 
zone the standard sampling Interval is one minute and the accumulation time Is 24 hours. 
Then the accumulation rises In discrete steps up to 4608 hours and the sampling interval 
Increases up to 48 minutes. The table Is organized with four orbits per page with the altitude 
constant and fixed Inclination orbits of 0°, 30°, 60°, and 90°. The first column under each 
inclination gives the dally fluence above the energy given In the first column of the page. The 
second column under each inclination gives the dally fluence In the energy band defined by the 
first two columns of the page. 
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Appendix A 
Function TRARA2 


Within Program MODEL and now MODEL87 the subroutine TRARA1 handles the interpolation 
of Z in E, L, and b, and so serves as the gateway to the numerical AE-8 models. A FORTRAN 
function, TRARA2, which is in TRARA1, handles the Z. L, and b interpolation linearly, but 
the paths are not along constants of the variables. TRARA2 has been discussed briefly in 
Teague et al. (1972), but this and the source listing for the function leave some points unclear 
about how TRARA2 really works. Because of this additional comments have been added to the 
source listing, and the names of some variables have been changed to make it easier to follow 
the logic. The original code was written in the days when the GOTO command was in vogue 
rather than structured programming. This has not been changed so that the function operates 
in exactly the same manner as before. This makes the code harder to follow, but the 
description of the method and the newly commented listing presented in this appendix should 
make everything clear to the user. 


Method 


The method can best be illustrated by using Figure Al. Here two Z versus x curves are shown 
for adjacent L values, LI and L2 (L2 > LI) at one of the model grid energies. Each grid 
point in Z-x space for LI and L2 is shown in Figure Al. The grid energy is fed to TRARA2 in 
the form of the starting address in the BLOCK DATA for that energy by the argument of MAP, 
which itself is an argument of TRARA2. Based on the input L value, TRARA2 determines LI 
and L2. The general idea is to And the two radial lines passing from the Z-x origin through 

any grid point that most closely brackets the input point. Since the Z coordinate of the input 
is what is desired, then an iterative approach is needed. The coordinates of the grid points are 
known: the equations of the two curves are known (these are straight lines connecting the grid 
points)' the equations of the radial lines are known. Except for OCB3 these radial lines are 
shown as line segments in Figure Al . The intersection of a radial line with the other curve can 
readily be computed iteratively. On these two radial lines the L interpolation is done to find 
the x values. Using this and the equation of the radial lines, the Z coordinates are obtained. 
Then the interpolation along the line defined by these two interpolated points provides the 
output Z value. 
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The specifics will be given by considering three Input points shown in Figure A1 as PI. P 2 , 
and P3. On the basis of the equatorial values, Z1 (LI ) and Z1 (L2 ) . TRARA2 sets up the' 
larger one as the top (T) set of 5x points and the other as the bottom (B) set. These sets are 
shown converted to x values in Figure A1 and are referred to as the T and B curves. Thus. 
Initially the LI curve is the T curve. The procedure to find the bracketing radial lines starts 
by stepping on the T curve from the equator to grid points until the grid point with an x value 
greater than the Input value FNx Is reached. 

Consider the PI case first. One can see that no grid point on LI satisfies this condition. Then 
the code reverses the role of LI and L2 so that L2 becomes the T curve. When the stepping 
reaches the grid point near PI labeled BT1, then the following procedure Is employed. The 
line OAT1 is constructed (shown as AB1 : ATI). Then the B curve (now LI) Is stepped from the 
equator until the grid point labeled CB1 Is found. The criterion used Is to compute the 
Intersection point AB1 between OAT1 and the B curve until 


x(ABl) < x(Bi) 

where i is the Index of the grid point. The formula for x (AB1) Is given by 


( A . 1 ) 


' Z ( AT 1 ) * (x (Bj) -x (Bj -t ) ) 
ZDEL*x (ATI ) 


+ 1 


(A. 2) 


Thus, the straight line that represents the B curve changes with the value of i. For the 
example being discussed, B 21 = CB1 and the equation of the line representing curve B are 
considerably different than for the segment AT 3 : BT3. One can see by inspection that any grid 
point on LI closer to the equator would not satisfy the condition. The assumption Is now 
made that the x values vary linearly In L along OAT1, so a linear interpolation in L Is done 
to obtain x(ATl). The Z(AT1) value is given by 


Z (All) 


Z (ATI) x (AB1 ) 
x (ATI) 


(A. 3) 


These coordinates are then stored in the XPCM and ZPCM registers, which are used to hold the 
trial coordinates for the nearest point with XPCM less than FNx. Next the slopes of OBT1 and 
OCB 1 are compared, and the larger one ( OCB 1 ) is used first to produce the point labeled C 1 1 
whose coordinates are placed In the XPCP and ZPCP registers. Then the check is made If 
x(CIl) Is greater than or equal to x(Pl). Since it Is not true In this case, the PCP registers 
are dumped Into the PCM registers (OCT1 replaces 0AT1 ). and the B curve is stepped and 
tested until the test slope is less than the OBT slope. In this case the next grid point tested 
would be that of BT1. Since x(BIl) Is greater than FNx, the final Interpolation In x can be 
done. This final Interpolation assumes that Z is linear in x. Knowing the coordinates of BI 1 
and CI1, the value of Z(P1) is computed. 

If one now considers P2, which has the same x value as PI. things would proceed In the same 
fashion as for PI up to the test x(BIl) Is greater than x(Pl). However, the position of P 2 
Implies L(P2) Is less than L(P1), so all the Interpolated points All, BI1, and CI1 would 
be closer to the origin than shown In Figure Al. Thus, the new BI1 would fall the test x(BIl) 
Is greater than x (P2 ) . This would result in testing the next grid point radial line that had the 
larger slope, and the ping-ponging between the LI and L2 curve would continue until the lines 
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0AT2 and 0BT2 would be reached. The final solution would be done In the same manner as 
described above. 

What has been described above Is the most straightforward way for the program to proceed. 
However, there is a conservative approach in the logic that produces a path that leads to the 
correct answer but apparently is not the most direct path. Tilings actually proceed as described 
above up to the place where ATI and BT1 are Identified. Because there has been a reversal of 
the curves, the logic now forces a path through the code that has been devised to handle the case 
where FNx is near the equator. This procedure makes OAT (the numeric suffix is dropped at 
this point) the Z axis so that the determination of A I does not involve any interpolation 
using x. Then BT is the first step on the T curve and CB is the first step on the B curve. The 
huger slope point is tested for the solution as described above, and the ping-pong procedure is 
started. In the case of FNx near 0, either this first step BT or CB provides the XPCP. In the 
P 1 example the failure of the first step to provide an XPCP starts the ping-pong procedure that 
results in the same solution as described previously. 

Now consider P 3 because it illustrates another aspect of the logic employed in the code. No 
reversal is needed because the point BT3 satisfies the condition. In this case the LI curve 
remains the T curve. Thus, AT 3 and BT3 are identified as shown in Figure Al. However, the 
fact that now x (AB 3) is greater than x (AT 3) causes the logic to force the computation to 
pass through the FNx near 0 procedure described above. Here again the ping-pong procedure 
produces the proper answer, which in this case uses the lines OAB3 and OBB3 shown near P3. 
If in Figure Al Z 1 ( L2 ) is greater than Z 1 ( L 1 ) , then for all cases considered above the 
procedure would progress as initially described above, except no reversal would be needed for 
airy of the cases considered. 

There are several other cases handled in the code, but these are quite transparent. It is believed 
the code is now commented well enough that users can understand in detail exactly what is 
being done. 


Source Code Listing 

FUNCTION TRARA2 (MAP, IL, IX) 

0 ***************************************************************** 


C*** TRARA2 INTERPOLATES LINEARLY IN L-B/BO-MAP TO OBTAIN *** 
C*** THE LOGARITHM OF INTEGRAL FLUX AT GIVEN L AND B/BO. *** 
C*** INPUTS: MAP ( . . ) IS SUB-MAP FOR SPECIFIC GRID ENERGY *** 
C*** OF TRAPPED RADIATION MODEL MAP *** 
C*** IL SCALED-L VALUE *** 
C*** IX SCALED- (B/BO - 1)=X VALUE *** 
C*** OUTPUTS: TRARA2 SCALED-LOGARITHM OF PARTICLE FLUX=Z *** 
C*** MOST WORKING VARIABLES, I.E., THOSE THAT CHANGE MEANING *** 
C*** IN COURSE OF THE PROGRAM ARE DENOTED WITH THE SUFFIX *** 
C*** LETTER, W. THIS DOES NOT APPLY TO INDEXES. *** 

Q************ *********************************************** ****** 

C*** SEE MAIN PROGRAM 'MODEL' OR ' RADBELT ' FOR EXPLANATION *** 
C*** OF MAP FORMAT SCALING FACTORS *** 
C*** THE PARAMETER THAT DEFINES THE FIXED CHANGE OF Z BE- *** 
C*** TWEEN X GRID POINTS IS OBTAINED FROM ' COMMON /TRA2/ • *** 

Q* *************************************************************** * 
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DIMENSION MAP (1 ) 

COMMON / TRA2 / F I STEP 
C 

C CHANGE INTEGER INPUTS TO FLOATING POINT, RENAME FISTEP, AND 
C INITIALIZE. 

C 

FNL=IL 

FNX=IX 

ZDEL=FISTEP 

ITIME=0 

IT=0 

C 

C FIND CONSECUTIVE SUB-SUB MAPS FOR SCALED-L VALUES, LSB & LST, 

C WITH LSB.LT.IL AND IL.LE.LST . LB & LT ARE LENGTHS OF THE 
C SUB-SUB MAPS. IB+1 & IT+1 ARE THE INDEX VALUES IN MAP 
C WHERE THE VALUES OF LB & LT ARE STORED. IB+2 & IT+2 ARE 
C THE INDEX VALUES WHERE LSB AND LST ARE STORED. IB+3 & IT+3 
C ARE THE INDEX VALUES WHERE THE EQUATORIAL Z VALUES, Z1B AND 
C Z1T ARE STORED. 

1 LT=MAP (IT+1 ) 

IF (MAP (IT+2) .GT. IL) GOTO 2 ! PROPER LST HAS BEEN FOUND 

IB=IT 
LB=LT 

IT=IT+LT ! CHECKING ON NEXT L VALUE 

GOTO 1 
2 CONTINUE 
C 

C IF A SUB-SUB MAP IS EMPTY, I.E., LENGTH LESS THAN 4, THEN 
C TRARA2=0; DOUBLE EMPTY SETS OCCUR AT THE INNER AND OUTER 
C BOUNDARIES OF THE MODELS. 

C 

IF( (LB.LT.4) .AND. (LT.LT.4) ) GOTO 50 ! THIS SETS TRARA2=0 . 

C 

C IF Z1T IS LESS THAN Z1B INTERCHANGE VARIABLES AND INDEXES SO 
C THAT THE SUFFIX LETTER T IS ASSOCIATED WITH THE L VALUE THAT 
C HAS THE LARGER Z1 . THIS MAKES IT EASIER TO FOLLOW THE LOGIC OF 
C THE CODE. 

C 

IF (MAP (IT+3) .GT. MAP (IB+3) ) GOTO 10 !NO INTERCHANGE NEEDED 
5 KT=IB 

IB=IT 
IT=KT 
KT=LB 
LB=LT 
LT=KT 


DETERMINE THE INTERPOLATION FACTOR IN SCALED-L VALUE 

10 FLLB=MAP (IB+2) 

FLLT=MAP (IT+2) 

DFL= (FNL-FLLB/ (FLLT-FLLB) 

Z1B=MAP (IB+3) 

Z1T=MAP (IT+3) 

XBW=0 . ! SET WORKING VARIABLES 

XTW=0 . 

ZBW=Z1B 

ZTW=Z1T 

IF (LB.LT.4) GOTO 32 ! CASE WHERE Z1B=0 . 
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c 

C THE X LOOPING IS BASED ON FINDING, FOR A GRID ENERGY AND THE 
C TWO SELECTED SCALED-L VALUES, IN X-Z SPACE TWO POINTS, PCM AND 
C PCP, WHERE X (PCM) AND X(PCP) BRACKET THE INPUT X VALUE, FNX. 

C EACH SUB-SUB MAP CONTAINS THE INFORMATION TO CONSTRUCT A PLOT 
C OF POINTS GIVING Z VERSUS X SO WILL BE REFERRED TO AS A CURVE. 

C FOR THE LINEAR INTERPOLATION USED THESE GRID POINTS ARE CONNECTED 
C BY STRAIGHT LINES. THERE ARE SEVERAL SPECIAL CASES BUT THE GENERAL 
C APPROACH IS TO START AT THE EQUATOR (X=0) FOR THE TOP CURVE AND 
C STEP FROM GRID POINT TO GRID POINT UNTIL FNX IS 
C BRACKETED. TAKING THE SMALLER XT GRID POINT, AT, ONE DRAWS A 
C RADIAL LINE, OAT, TO THE X-Z ORIGIN. THE INTERSECTION OF OAT 
C AND THE LOWER CURVE IS POINT AB. IT CAN BE DETERMINED DIRECTLY 
C SINCE THE EQUATIONS OF BOTH LINES ARE KNOWN. THE SECOND XT GRID 
C POINT IS BT AND THE INTERSECTION OF LINE OBT WITH B CURVE IS BB. 

C A GRID POINT ON THE B CURVE IS CB AND THE INTERSCTION OF LINE OCB 
C WITH THE CURVE IS CT. USING COMPUTED X (AB) AND X (CT) 

C AN L-INTERPOLATED X (PCM) IS FOUND. THEN ONE ITERATES BY STEPPING 
C IN X FIRST ON THE B CURVE AND THEN IF NEEDED ON THE T CURVE IN A 
C PING PONG FASHION TO GET THE BRACKETING. THEN AN X, Z 
C INTERPOLATION USING PCM AND PCP IS MADE. FINALLY THE INTERPOLATED 
C Z VALUE IS RETURNED. 

C 

DO 17 J2=4,LT 

XT I NCR=MAP (IT+J2) 

IF (XTW+XTINCR.GT.FNX) GOTO 23 ! USUAL CASE, SELECTED POINT, AT 
XTW=XTW+XTINCR 
17 ZTW=ZTW-ZDEL 

ITIME=ITIME+1 ! SETS FLAG THAT NO SUCH POINT FOUND ON CURVE. 

C 

C IN THIS CASE FNX LIES IN THE RANGE BETWEEN CUTOFFS OF THE TWO 
C L VALUES AND REVERSAL OF CURVES IS REQUIRED. 

C 

IF (ITIME.EQ. 1) GOTO 5 ! LEADS TO REVERSAL OF T AND B CURVES 
GOTO 50 ! SET TRARA2=0 . REVERSAL FAILED. 

23 IF (ITIME.EQ. 1) GOTO 30 
C 

C THE GOTO 30 OCCURS FOR THE CASE WHERE THERE IS A REVERSAL AND THE 
C POINT, AT, HAS BEEN FOUND. IT IS NOT CLEAR WHY THE CODE IS WRITTEN 
C THE WAY IT IS. PROCEEDING AT 30 RESETS THE LINE OAT TO BE THE Z 
C AXIS (X=0 ) AND THINGS PROCEED LIKE FNX IS NEAR THE EQUATOR. THE 
C TEST IN 15 THROWS THIS OUT AS A SOLUTION AND THE PING PONG BETWEEN 
C 15 AND 20 EXPLAINED BELOW RESULT IN STEPPING OUT ON BOTH CURVES 
C UNTIL THE PROPER TWO BRACKETING LINES ARE FOUND. THUS THE RIGHT 
C ANSWER IS OBTAINED BUT THE LOGIC IS CONFUSING. 

C 

IF ( J2 . EQ . 4 ) GOTO 28 !THIS IS CASE WHERE FNX IS NEAR 0. 
SLTW=ZTW/XTW ! SLOPE OF RADIAL LINE, OAT 
C 

C THE GRID POINT, AT, HAS BEEN FOUND. NOW THE B CURVE IS STEPPED 
C AND THE X (AB) =XABW BELOW IS COMPUTED UNTIL ITS VALUE IS LESS THAN 
C THE VALUE OF X(CB)=XBW BELOW. 

C 

DO 27 Jl=4, LB 

XBINCR=MAP (IB+J1) 

XBW=XBW+XBINCR 

ZBW=ZBW-ZDEL 


A-6 



Appendix A 


XABW= ( (ZBW/ZDEL) *XBINCR+XBW) / ( (XBINCR/ZDEL) *SLTW+1 . ) 

IF (XABW.LE.XBW) GOTO 31 
27 CONTINUE 

IF (XABW.LE.XTW) GOTO 50 ! TRARA2=0 . 

31 IF (XABW.LE.XTW) GOTO 29 ! USUAL CASE . ONLY FALSE FOR REVERSAL 
C 

C 31 WILL NOT BE REACHED IF THERE IS A REVERSAL OR FNX IS NEAR 0. 

C FAILURE OF TEST IN 31 OCCURS ONLY WHEN T CURVE LIES BELOW B CURVE 
C IN THE REGION OF THE INPUT VALUE FNX. IT IS NOT CLEAR 
C WHY SUCH CASES ARE ROUTED THIS WAY. HOWEVER, THE 

C CORRECT ANSWER IS OBTAINED. CODE BELOW TREATS PROPERLY THE CASE 
C WHERE FNX IS NEAR 0. THUS OAT = Z AXIS AND ZPCM IS COMPUTED USING 
C ONLY L INTERPOLATION. THE SLOPES OF OBT AND OCB ARE EXAMINED TO 
C SELECT THE LARGER ONE AT 15. THEN XPCP AND ZPCP ARE COMPUTED IN 15 
C OR 20. IF THIS IS A SOLUTION, I.E. XPCP. GE. FNX, THE FINAL ANSWER 
C IS COMPUTED. IF NOT IT STEPS THE CURVE AND TESTS UNTIL THE TEST 
C SLOPE IS LESS THAN THE OTHER CURVE'S SLOPE. 

C 

XBW=0 . ! SETTING TO EQUATOR 

30 XTW=0 . ! AT FOUND AFTER A REVERSAL 

32 J2=4 ! CASE FOR Z1B=0 

XTINCR=MAP (IT+J2) 

ZBW=Z1B 

ZTW=Z1T 

28 ZPCM=Z1B+ (ZTW-ZBW) *DFL !Z VALUE AT X=0;CASE FOR FNX NEAR 0 
XPCM=0 . ! X VALUE FOR ZPCM 

XTW=XTW+XTINCR ! 1ST STEP ON T CURVE, BUT IS BT NOW 

ZTW=ZTW-ZDEL 

SLTW=ZTW/XTW ! OBT SLOPE 

IF (LB. LT . 4 ) GOTO 35 ! LB IS AN EMPTY L SET SO XBINCR=0 . 

Jl=4 

XB I NCR=MAP (IB+J1 ) 

XBW=XBW+XBINCR ! 1ST STEP ON B CURVE, CB 

ZBW=ZBW-ZDEL 
S LBW= Z BW / XBW 
GOTO 15 
C 

C THIS GENERATES PARAMETERS FOR THE NON-NEAR EQUATORIAL CASES 
C 

29 XPCM=XABW+ (XTW-XABW) *DFL !L INTERPOLATION BETWEEN AT & AB 
ZPCM=XPCM*SLTW 

ZTW=ZTW-ZDEL !NOW THE T CURVE IS STEPPED. 

XTW=XTW+XTINCR 
SLTW=ZTW/XTW ! SLOPE OF OBT 
SLBW=ZBW/XBW ! SLOPE OF OCB 
C 

C THE CODE IN 20 COMPUTES INTERSECTION POINTS ON THE B CURVE AND 
C STEPS ALONG THE T CURVE UNTIL THE SLOPE OF THE RADIAL LINE TO THE 
C GRID POINT IS LESS THAN THE SLOPE OF THE RADIAL LINE THROUGH THE 
C LAST SELECTED POINT ON THE B CURVE. THE CODE IN 15 HANDLES THE 
C SAME FUNCTIONS WITH THE B AND T CURVES REVERSED. IF THE PATH TO 
C 60 IS FROM 15 CODE THEN XPCP COMES FROM A T GRID RADIAL LINE. A 
C B GRID RADIAL LINE COMES FROM 20 CODE. THE SELECTION AT 15 ALWAYS 
C MAKES THE LARGER SLOPE RADIAL LINE THE FIRST TO BE TESTED. IF 
C THAT FAILS THAT CURVE IS STEPPED AND TESTED UNTIL THE PCP POINT 
C IS FOUND OR THE TEST SLOPE BECOMES SMALLER THAN THE SLOPE OF THE 
C LINE THROUGH THE GRID POINT WAITING TO BE TESTED ON THE OTHER 
C CURVE. THIS IS THE PING PONG. 
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15 IF(SLBW.LT.SLTW) GOTO 20 

XCTW= ( (ZTW/ZDEL) *XTINCR+XTW) / ( (XTINCR/ZDEL) *SLBW+1 . ) 
XPCP=XBW+ (XCTW-XBW) *DFL ! L INTERPOLATION 
ZPCP=XPCP*SLBW 

IF (XPCP.GE.FNX) GOTO 60 .'READY FOR FINAL X INTERPOLATION 
C 

C MUST BEGIN ITERATION TO GET POINTS PCM AND PCP BRACKETING INPUT 
C SO PCP REPLACES PCM AND NEXT POINT ON B CURVE IS EXAMINED TO FIND 
C A PCP. 

C 

XPCM-XPCP 
ZPCM-ZPCP 

IF(Jl.GE.LB) GOTO 50 ! LB GRID IS EXHAUSTED, SET TRARA2=0 . 
J1=J1+1 

XB I NCR=MAP (LB+ J1 ) 

ZBW=ZBW-ZDEL 
XBW=XBW+XBINCR 
SLBW=ZBW/XBW 

GOTO 15 ! NOW RETURN TO SEE IF 20 OR 60 EXIT IS OPEN 

THE BB OR SUBSEQUENT INTERSECTION POINTS ARE COMPUTED AND TESTED 

20 XBBW- ( (ZBW/ZDEL) *XBINCR*+XBW) / ( (XBINCR/ZDEL) *SLTW+1 . ) 
XPCP=XBBW+ (XTW-XBBW) *DFL 
ZPCP-XPCP*SLTW 

IF (XPCP.GE.FNX) GOTO 60 ! USUAL CASE 
C 

C NOW THE TOP CURVE MUST BE STEPPED TO THE NEXT POINT AND THE 
C PCP VALUES MUST BE TRANSFERRED TO THE PCM CELLS 
C 

XPCM=XPCP 

ZPCM=ZPCP 

IF ( J2 .GE . LT) GOTO 50 ! END OF TOP CURVE, SET TRARA2=0 . 

J2=J2+1 

XTINCR=MAP (I2+J2) 

ZTW=ZTW-ZDEL 
XTW=XTW-XTINCR 
SLTW=ZTW/XTW 
GOTO 15 
35 XBINCR=0 . 

SLBW=-900000. ! T CURVE WILL ALWAYS BE TESTED AT 15 BRANCH 

GOTO 20 ! ENTERS 20 WITH XBW=0 . 

60 IF (XPCP . LT . XPCM+1 . E-20 ) GOTO 50 

TRATRA2=ZPCM+ (ZPCP-ZPCM) * ( (FNX-XPCM) / (XPCP-XPCM) ) 

TRARA2 =AMAX 1 (TRARA2, 0 . ) 

RETURN 
50 TRARA2=0. 

RETURN 

END 
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Table 1. Azur Greater Than 1.5 MeV Normalized Local Time Distribution 


Local Time 

L = 3 

L = 4 

in 

II 

A 

L = 5.5 

L = 6 

(Hrs.) 

<D/K 

<D/K 

O/K 

O/K 

O/K 


1 


2 

— 

— 

0.767 

0.672 

0.698 

3 

— 

0.700 

0.894 

0.810 

0.594 

4 

1.038 

0.795 

0.853 

0.992 

1.081 

5 

1.122 

0.849 

1.024 

0.960 

0.861 

6 

0.968 

0.973 

0.951 

0.885 

0.796 

7 

0.953 

0.996 

1.026 

1.034 

0.951 

8 

0.820 

1.039 

1.087 

1.064 

1.257 

9 

— 

1.130 

1.127 

1.254 

1.203 

10 

- 

- 

- 

1.151 

1.385 

11 

- 

- 

- 

- 

- 

12 

— 

- 

- 

- 

- 

13 

— 

— 

- 

1.024 

1.434 

14 

— 

— 

1.099 

1.160 

1.170 

15 

— 

1.055 

1.057 

1.065 

1.199 

16 

0.952 

1.033 

1.096 

1.231 

1.048 

17 

0.929 

0.992 

1.051 

1.045 

1.040 

18 

1.079 

1.057 

1.047 

0.987 

0.961 

19 

1.014 

1.063 

1.075 

1.071 

1.031 

20 

0.968 

0.997 

1.008 

0.940 

1.002 

21 

— 

0.889 

0.811 

0.904 

0.916 

22 

— 

- 

— 

0.629 

0.859 

23 

- 

- 

- 

- 

- 

24 

- 

- 

- 

- 

- 


C 

- 0.049 

0.053 

0.063 

0.090 

0.106 

95% Confidence 
Interval Half 

0.003 

0.004 

0.003 

0.005 

0.007 

Width 

K 

0.981 

0.962 

0.992 

0.970 

0.991 
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Table 3. Ratio of (CR-BREM)/CR for OV1-19 Spectrometers 


Detector 

ID 

(keV) 




T. 





3 

3.5 

4 

4.5 

5 

5.5 

6 

6.6 

LEI 

53 

-0.83 

-3.56 

-4.42 

-1.86 

-1.30 

-0.63 

-0.39 

-0.75 

LE2 

52 

0.92 

-0.55 

-0.51 

0.63 

0.69 

0.83 

0.88 

0.80 

LE4 

192 

0.99 

0.87 

0.92 

0.98 

0.99 

1.00 

1.00 

0.99 

LE5 

250 

0.99 

0.90 

0.91 

0.97 

0.99 

0.99 

0.99 

0.99 

LE6 

312 

0.97 

0.79 

0.77 

0.93 

0.96 

0.98 

0.98 

0.97 

LE7 

376 

0.95 

0.71 

0.67 

0.89 

0.95 

0.97 

0.97 

0.95 

LE8 

444 

0.92 

0.56 

0.53 

0.85 

0.92 

0.95 

0.95 

0.90 

HEl 

537 

0.97 

0.96 

0.98 

0.98 

0.98 

0.99 

0.97 

0.98 

HE2 

822 

0.51 

0.57 

0.76 

0.83 

0.86 

0.90 

0.85 

0.41 

HE5 

1711 

-2.12 

0.79 

0.80 

0.61 

0.34 

0.06 

-1.09 

-0.53 

HE6 

2018 

-0.86 

0.73 

0.70 

0.48 

0.14 

-0.06 

-1.54 

-0.24 

HE7 

2321 

-6.62 

0.27 

0.10 

-1.29 

-3.75 

-6.24 

-14.03 

-3.69 

HE8 

2632 

-2.69 

0.50 

0.24 

-1.03 

-2.80 

-4.53 

-7.72 

-1.90 

HE9 

2932 

-7.55 

0.16 

-0.70 

-3.55 

-7.49 

-9.49 

-7.50 

-3.06 

HEl 1 

3557 

-8.09 

-1.76 

-4.40 

-11.40 

-14.79 

-12.77 

-6.22 

-2.20 

HE12 

3863 

-9.82 

-1.12 

-3.56 

-9.15 

-14.67 

-16.70 

-9.35 

-3.76 

HE13 

4169 

-13.97 

-3.08 

-8.00 

-17.57 

-22.13 

-20.16 

-6.87 

-3.53 

HE14 

4476 

-9.13 

-2.04 

-5.15 

-10.62 

-13.60 

-12.51 

-5.68 

2.46 

HE15 

4783 

-11.46 

-3.61 

-8.63 

-16.49 

-19.33 

-16.22 

-6.42 

-3.16 

HE16 

5091 

-12.36 

-3.33 

-8.80 

-16.30 

-18.61 

-18.77 

-7.41 

-3.18 


Table 4. Flux Fraction d for OV1-19 HES Channels 


L d Remarks 


3.0 

0.1450 

Determined using AZUR 

3.2 

0.1040 

Determined using AZUR 

3.4 

0.0744 

Determined using AZUR 

3.5 

0.0720 

Determined using AZUR 

3.6 

0.0549 

Determined using AZUR 

3.8 

0.0444 

Determined using AZUR 

4.0 

0.0486 

Determined using AZUR 

4.2 

0.0496 

Determined using AZUR 

4.4 

0.0503 

Determined using AZUR 

4.5 

0.0440 

Determined using AZUR 

4.6 

0.0541 

Determined using AZUR 

4.8 

0.0551 

Determined using AZUR 

5.0 

0.0547 

Determined using AZUR 

5.5 

0.0478 

Determined using AZUR 


> 4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

4500 keV data. 

L = 5.5 and 6.6; d AZUR = .0527. 


6.0 0.0322 Determined by interpolation between 

6.6 0.0145 Determined using ATS 6 > 3900 keV data; d AZUR = .0592. 


> 

> 

> 

> 

> 

> 

> 

> 

> 

> 

> 

> 

> 
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Table 5. Parameters Cor Conversion of OV1-19 Flux to Integral Energy Spectrum 


Electron Energy (keV) 

Useful Channel 


Channel 

Index-k 

Detector 

ID 

Upper 

Limit 

Effective 

Midpoint 

Lower 

Limit 

Width 

W(k) 

1 

LEI 

71 

55 

39 

32 

2 

LE2 

140 

106 

71 

69 

3 

LE4 

220 

180 

140 

80 

4 

LE5 

281 

251 

220 

61 

5 

LE6 

344 

313 

281 

63 

6 

LE7 

409 

377 

344 

65 

7 

LE8 

448 

429 

409 

39 

8 

HE1 

677 

563 

448 

229 

9 

HE2 

1266 

972 

677 

589 

10 

HE5 

1866 

1566 

1266 

600 

11 

HE6 

2172 

2019 

1866 

306 

12 

HE7 

2478 

2325 

2172 

306 

13 

HE8 

2785 

2632 

2478 

307 

14 

HE9 

3237 

3011 

2785 

452 

15 

HE11 

3708 

3473 

3237 

471 

16 

HE12 

4016 

3862 

3708 

308 

17 

HE 13 

4325 

4171 

4016 

309 





Table 6A. OV1-19 Equatorial Unidirectional Flux Uncorrected for Background (electrons/ sq cm-ster-keV-s) 
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Table 8. Grid Points and Other Parameters for AE-8 Numerical Representation 


Energy Value 
(MeV) 

L Value 

L Value 

Other Values of Other 

Parameters Parameters 

Notes/Remarks 

0.04 

0.00 

3.20 

Model Names 

AE-8MIN 


0.10 

1.10 

3.40 


AE-8 MAX 


0.25 

1.20 

3.50 

Total Storage 

13168 

Cells for AE-8MIN. 

0.50 

1.25 

3.60 


13548 

Cells for AE-8MAX. 

0.75 

1.30 

3.80 

Scale Factors 


Stored Value/Scale Factor = True Value. 

1.00 

1.35 

4.00 

Energy 

6400 


1.50 

1.40 

4.20 

L 

2100 


2.00 

1.45 

4.40 

b Increments 

1024 

Occurs when Z has dropped by 0.25. 

2.50 

1.50 

4.50 

Z1 

1024 

This parameter ZDEL = 0.25 does not 

3.00 

1.60 

4.60 



appear in BLOCK DATA statement but 

3.50 

1.70 

4.80 



is in Program MODEL87 that handles 

4.00 

1.80 

5.00 



BLOCK DATA statement. 

4.50 

1.90 

5.50 




5.00 

2.00 

6.00 




5.50 

2.20 

6.60 




6.00 

2.40 

7.00 




6.50 

2.50 

7.50 




7.00 

2.60 

8.00 





2.69 

- 




- 

2.70 

9.00 





2.79 

- 




- 

2.80 

10.00 





2.89 

- 




- 

2.90 

11.00 




- 

3.00 

12.00 




- 

- 

15.60 
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Ta ble 9. AE-8MAX Header and Logical Data Record (partial) 


C*«**»» AE8MAX 

BLOCK DATA 

DOUBLE PRECISION NAME 

COMMON/AE8MAX /NAME,EPOCH,SCAE,SCAL,SCAB,SCAF .LENGTH 
X.MAP 1(163), MAP 2(1 63), MAP 3(163),MAP 4(163),MAP 5(163) 

X.MAP 6(163), MAP 7(163),MAP 8(163),MAP 9(163),MAP 10(163) 

X MAP 1 1(163), MAP 12(163), MAP 13(163),MAP 14(163),MAP 15(163) 

X MAP 16(163), MAP 17(163),MAP 18(163),MAP 19(163),MAP 20(163) 

x’mAP 21(163),MAP 22(163),MAP 23(163),MAP 24(163),MAP 25(163) 

x!mAP 26(163) .MAP 27(163),MAP 28(163),MAP 29(163),MAP 30(163) 

X MAP 31(163), MAP 32(163),MAP 33(163),MAP 34(163),MAP 35(163) 

x’MAP 36 ( 1 63), MAP 37(163),MAP 38(163),MAP 39(163),MAP 40(163) 

xiMAP 41(163), MAP 42 (163), MAP 43(163),MAP 44(163),MAP 45(163) 

X.MAP 46(163), MAP 47(163), MAP 48(163),MAP 49(163),MAP 50(163) 

xiMAP 51(163), MAP 52(163), MAP 53(163),MAP 54(163),MAP 55(163) 

x!mAP 56(163),MAP 57(163), MAP 58(163), MAP 59(163),MAP 60(163) 

X MAP 61(163),MAP 62(163).MAP 63(163), MAP 64(163),MAP 65(163) 

X MAP 66(1 63), MAP 67(163),MAP 68(163).MAP 69(163),MAP 70(163) 

X MAP 71(163), MAP 72(163),MAP 73(163),MAP 74(163),MAP 75(163) 

X’MAP 76(1 63), MAP 77(163),MAP 78(163),MAP 79(163),MAP 80(163) 

xiMAP 81 (163), MAP 82(163), MAP 83(163),MAP 84(19) 

DATA NAME /8HAE8-MAX / 


DATA EPOCH / 

1990.000/ 






DATA SCAE / 

6400.000/ 






DATA SCAL / 

2100.000/ 






DATA SCAB / 

1024.000/ 






DATA SCAF / 

1024.000/ 






DATA LENGTH/ 


13548/ 






DATA MAP 1/ 








X 

1442, 

256, 

3. 

0. 

0 , 

3. 

2310, 

0 , 

27, 

X 

2520, 

5907, 

37, 

19, 

19, 

19. 

18, 

17, 

16, 

X 

15, 

13, 

12, 

11. 

10. 

9. 

8, 

7. 

6, 

X 

6, 

5, 

4, 

4, 

3, 

3, 

3. 

2, 

32, 

X 

2625, 

7393, 

72. 

33, 

33, 

32, 

30, 

28. 

26, 

X 

24. 

22. 

19, 

17, 

16, 

14. 

12. 

11, 

9, 

X 

8. 

7, 

6, 

6, 

5, 

4, 

4, 

3. 

2, 

X 

3. 

2, 

2. 

2. 

34, 

2730, 

7926, 

171. 

55. 

X 

49. 

46. 

42, 

38. 

34, 

30, 

27. 

23, 

20, 

X 

18, 

15. 

13, 

11. 

10, 

8. 

7. 

6, 

5, 

X 

4, 

4, 

2, 

3. 

2, 

1. 

2. 

2. 

1, 

X 

1 . 

0, 

35. 

2835, 

8153, 

279, 

85, 

72. 

65, 

X 

58. 

51, 

44, 

37. 

32, 

26, 

22, 

18, 

15. 

X 

12. 

10. 

8, 

6, 

6, 

4, 

4. 

2, 

3. 

X 

3, 

2. 

1 . 

0 , 

0, 

0, 

0. 

0. 

0 . 

X 

0. 

36. 

2940, 

8339, 

389, 

123, 

105, 

93. 

80, 

X 

67, 

55. 

45, 

36, 

29, 

23, 

18, 

14. 

11. 

X 

9. 

7. 

5, 

4, 

4. 

2, 

1 . 

1 . 

0. 

Y 

0/ 










00000010 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

00000080 

00000090 

00000100 

00000110 

00000120 

00000130 

00000140 

00000150 

00000160 

00000170 

00000180 

00000190 

00000200 

00000210 

00000220 

00000230 

00000240 

00000250 

00000260 

00000270 

00000280 

00000290 

00000300 

00000310 

00000320 

00000330 

00000340 

00000350 

00000360 

00000370 

00000380 

00000390 

00000400 

00000410 

00000420 

00000430 

00000440 

00000450 

00000460 

00000470 

00000480 
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Table 10. AE-8 Local Time Amplitude Function 
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Table 12. AE-8 Standard Deviation of the Z Function, 
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Table 13A. AE-8MAX Equatorial Omnidirectional Fluxes Above Energy E 
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Table 13B. AE-8MPT Equatorial Omnidirectional Fluxes Above Energy E 

V 0 SO MeV 0.75 MeV l.OOMeV l.SOMeV 2.00 MeV 2.50MeV 3.00 MeV 3.50 MeV 4.00 MeV 4.S0MeV 5.00 MeV 5.50 MeV 
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Table 14. Four Circular Orbit Dally Fluence S ummar y Using AE-8MAX 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Four Circular Orbit Daily Fluence Summary Using AE-8MAX (cont.) 
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Table 14. Pour Circular Orbit Dally Fluence Summary Using AE-8MAX (cont.) 
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Figures 
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Figure 1. Local Time Variations for B > 1.5 MeV 
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Magnetic Cutoff Parameters 
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Figure 11. 0V1-19 Data at L = 6.0 
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Equatorial Omnidirectional Flux > E (e/sqcm-a) 
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Figure 18. Comparison of AE-8MAX with Outer Zone Data at L = 3.0 
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Figure 19. Comparison of AE-8MAX with Outer Zone Data at L = 3.2 
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Figure 20. Comparison of AE-8MAX with Outer Zone Data at L = 3.4 
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Figure 23. Comparison of AB^IIAX with Outer Zone Data at L>&8 
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Figure 25. Comparison of AE-8MAX with Outer Zone Data at L = 4.2 
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Figure 26. Comparison of AE-8MAX with Outer Zone Data at L = 
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Figure 27. Comparison of AE-8MAX with Outer Zone Data at L = 
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Figure 28. Comparison of AE-8MAX with Outer Zone Data at L = 
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Equatorial Omnidirectional Flux > E (e/sqcm-*) 
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Figure 29. Comparison of AE-8MAX with Outer Zone Data at L = 4.8 
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Figure 32. Comparison of AB-8 with Outer Zone Data at L * 6.0 
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Figure 33a. Comparison of AE-8 with Some Outer Zone Data at L = 6.6 
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Figure 33b. Comparison of AE-8 with Some Outer Zone Data at L = 6.6 
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Figure 35. Comparison of AE-8 and AE-4 at L = 7.5 
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Figure 36. Comparison of AE-8 and AE-4 at L * 8.0 


F-36 




Equatorial Omnidirectional Flux > E (e/aqem-a) 










Equatorial Omnidirectional Flux > B(e/*q cm-*) 




Equatorial O mnidir ectional Flux > E(e/sq cm-*) 




Equatorial Omnidirectional Flux > E(e/tq cm-a) 


□ 0.25Max 



O 2. 00 Max 


A 4.50Max 


□ 6.50Max 


■ 0.25Min 


• 2.00Min 


• 4.50Min 


♦ 6.50Min 


Figure 43 . AE-8MAX/MIN Comparison at 0 . 25 , 2 . 00 , 4 . 50 , and 6.50 MeV 
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Figure 45. AE-8MAX/MIN Comparison at 0.75 and 3.00 MeV 
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